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(57) Abstract: The invention is directed to a method for diagnosing the state of a system. The system may be mechanical, chemical, 
electrical, medical, industrial, business operations, manufacturing related, and/or processing related, among other. The method may 
measure a signal from the system. Further, the method may compare the signal to an expected signal. The method may calculate a 
signal strength and/or a noise. The signal strength and noise may be functions of a frequency. Further, the signal strength and noise 
may be used to determine a channel capacity and/or a rate of information. A comparison ot the rate of information and the channel 
capacity may yield information associated with the state of the system. The information may be used in diagnosing the state of the 
system. Further, the expected signal may be derived from a model. The model may be tuned to the measured signal. The model may 
have parameters that are associated with features and/or faults of the system. These parameters may be used in diagnosing the state 
of the system. Further, selectively repeated diagnosis over time may yield a prognosis of the system. 
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MODEL-BASED MACHINE DIAGNOSTICS AND PROGNOSTICS USING THEORY OF 
NOISE AND COMMUNICATIONS 

5 BACKGROUND OF THE INVENTION 

Related Applications 

This application claims priority of U.S. patent Application, Serial No.60/235,251, filed 
September 25, 2001 entitled: "MODEL-BASED MACHINE DIAGNOSTICS AND PROGNOSTICS 
USING THEORY OF NOISE AND COMMUNICATIONS", and is incorporated herein by reference in 
10 its entirety. 

1. Field of the Invention 

The present invention generally relates to a method of diagnosing systems. In particular, the 
present invention relates to a model-based method tor diagnosing the operational health of a system, 
15 and for forecasting the future operational health of the system. 

2. Description of Prior Art 

Diagnostics and prognostics are used in many fields. These fields may include mechanical, 
chemical, electrical, medical, manufacturing, processing, and business operations, among others. Each 
of these fields has problems and difficulties relating to determining the source of a problem, identifying 
20 the severity of the problem, and predicting the behavior of a system in relation to the problem 

For example, reliability and maintenance of complex equipment is critical to productivity and 
product quality. The purchase price of many typical equipments may account for half the equipments' 
cost Maintenance and support during the "life-cycle" may consume an amount roughly equal to the 
book value of the asset 

25 For example, at a typical chip plant, billions of dollars are invested in equipment many 

traditional manufacturing plants invest hundreds of millions of dollars in manufac turing rnachinery. 
Most maintenance is rigidly scheduled by "time in service", not condition. Machine productivity is lost 
during maintenance downtime. This and unscheduled downtime due to failures represent a very large 
part of a rnachine user's cost of operation. 

30 In these ernbodiments, machines arc complex systems of components: gears, shafts, bearings, 

motors, lead screws, sensors, electronics, microprocessors, etc. integrated into a working whole. 
Machines in this context can also be biological, chemical, or hydraulic, among others. Defective or 
degraded components, alone or interacting, can render a machine dysfunctional. The rnachine may fail 
catastrophically and not complete its task, or it may lose tolerance, resulting in defective parts. 

35 Although methods and models exist for many individual component failures, errors from shghtly 
degraded components can "stack", yielding overall system failure even when these models predict 
health of all individual components. 
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Many typical designers of machinery and engineers that maintain machinery, essentially know 

what the potential system faults are, and at what locations in the machine these faults will occur. 
However, unexpected breakdowns will still occur. Many typical designers and engineers do not know 
and/or can only poorly predict when faults will occur. Further, periodically, healthy machinery must be 
5 taken out of service for maintenance. Thus adding an unnecessary cost Also, it may be very difficult 
and sometimes impossible to observe many of the conditions internal to the machine that lead to failure. 

One can easily imagine rnetaphoric extensions of these problems to other fields such as 
chemical, electrical, medical, manufacturing, processing, and business operations, among others. As 
such, many typical systems suffer from deficiencies in providing accurate diagnostics and prognostics. 
10 Many other problems and disadvantages of the prior art will become apparent to one skilled in the an 
after comparing such prior an with the present invention as described herein. 
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SUMMARY OF THE INVENTION 

Aspects of the invention are found in a method to diagnose operational health of a system, and 
to forecast fixture health. For example, the method may permit intelligent scheduling of maintenance 
downtime in a mechanical or chemical system. Further, the method may be used, for example, to avoid 
5 functional and catastrophic failures. 

Further aspect of the invention may be found in the method assembling models of the machine 
system, including system components and known system faults. Faults may be treated as "noise". In 
addition, parameters in the model may be "tuned" from signals from the real system, causing the model 
to mirnic the real system in its present condition Diagnosis may then be performed by observing the 
10 model. 

In another aspect of the invention, the method may treat the system as a communications 
channel, estimate signal and noise levels, and diagnose health of the system with a tuned model by 
assessing how much information per unit time the system in its present condition can convey over its 
communications channel. The method may compare this maximum amount to the amount required by 
15 the system to execute a certain task. The method may assess if a system, in a given state, can perform a 
certain function, with a specified performance, within a given tolerance. 

Other aspects of the invention may be found in a method based on fundamental principles of 
physics and urformanon theory. Further aspects may be found in the method assessing functional 
condition or state operable to perform a specified task, in addition to potential for catastrophic failure. 
20 Additionally, the method may operate on a tuned model, avoiding interpretation of complicated signals. 
Furthermore, the method may allow predictive scenarios for a system's possible future health and 
functional condition, given certain observed trends. 

In another aspect, for a different system or a new system design, only the model may be altered 
and not the basic diagnostic algorithm The method may also permit the incorporation of knowledge of 
25 faults, and the intent of the designers of systems, into the diagnostics routines. 

Another aspect of the invention maybe found in assembling detailed dynamic systems models 
of the system in question. The models may posses a one to one correspondence between at least a 
portion of the physical components or elements in the real physical system, and elements in the dynamic 
systems model. In some embodiments, all of the physical components are modeled in a one to one 
30 correspondence. 

In one embodiment the model may include all possible faults and potential failures in the 
system models. These effects may be tabulated as "noise" in the system Noise in a signal is the 
difference between the actual signal and the expected signal. In the model, noise may be induced by 
changes in parameters of dynamic system elements, which then alters any signals passing through a 
35 system. Alternatively, if a certain fault cannot be described by these means, then sources of noise may 
be inserted into the system model, at locations in the model that are consistent with the locations of the 
faults in the real system. The intensity of these noise sources may be adjusted to make the model behave 
like the real system 
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In an exemplary embodiment, the method may further include, placing sensors on a machine, 
to monitor the machine; exciting the machine; and observing the machine's response via the sensor 
outputs. The collected data may be used to tune the model's parameters, so the model immics the real 
system After data has been collected on the actual machine, the system model may be excited with the 
same excitation as the actual machine. The outputs of the model may be compared to the corresponding 
outputs of the real machine. If the model's outputs differ from the real machine's outputs, me values of 
model parameters may be adjusted or changed, including the intensity of the noise sources, until the 
model's outputs approximate the actual system's outputs. 

In one embodiment, the channel capacity, C, of the system may be calculated The channel 
capacity may be the maximum amount of information per unit time that can be measured from 
successfully conveyed through the machine. The channel capacity may depend on the design and 
construction of the system, and the present condition of the system, which results from aspects. These 
aspects may include manufacturing, aging and damage, among others. For example, faults may be 
encoded as "noise" in the modeL Analytically, the channel capacity may depend on the strength of the 
15 noise levels in the system relative to the strength of the excitation system response signal. 

In the exemplary embodiment, for a desired task to be performed by the machine, the rate of 
infonraiion associated with the task may be calculated The rate-of information may depend on the 
desired speed at which the machine does the task, the desired loads imposed on the machine, the 
complexity of the task, and the desired accuracy at which the machine should do the task Further, the 

20 rate of information may be measured 

Another aspect of the invention may be found in comparison of the rate of irtformation to the 
channel capacity. This comparison may be used to evaluate the operabiHty of the system If the rate of 
ir.formarion is less than or equal to the channel capacity, the system may perform the desired task within 
the desired precision. If the rate of Mormation is greater than the channel capacity, the system may 

25 functionally fail. 

Another aspect of the invention may be found in the loimulation of extremely detailed models 
of the system to describe a system's condition. In one exemplary emtotiurient of a system, the model 
includes bond graph based models of a motor, a gear box, and other mechanical transmission 
components. These extremely detailed models (a) exhibit a one to one correspondence between 

30 elements in the model and components in the real system; (b) incorporate many typical effects of the 
device into the model, including defects; (c) include in the models, via finite element concepts instilled 
into bond graphs, the dynamically distributed nature of components in the real system, and (d) use noise 
sources to account for defects and degradation of components. Simulation of the motor and gear box 
models may generate the complex spectra measured during operation of these devices. These models 

35 may irnraic real system behavior and may be used to store information regarding the health condition of 
the machine. 
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In a further aspect, the models tabulate the effects of system faults (system maladies) as 
"noise" in the machine. Noise may be the difference between the actual signal received, and the 
expected signal that should be received. As a machine degrades or ages, the difference between actual 
and expected signals may become larger. Thus noise levels may increase. These noise methods permit 
5 incorporation of faults into me models that heretofore could not be described analytically. The herein 
described methods have imported this body of knowledge to mechanical, hydraulic, other physical 
systems, and others, to name a few. 

In an additionally aspect, the method may be used to predict the future conditions of systems, 
for scheduling maintenance and avoiding functional and catastrophic failures of the systems. The 
10 method may forecast if a complex system is capable of performing a given task, at a given speed and 
load, within a specified tolerance. 

The model system may be implemented on a computer system. Hardware and software 
components may in combination allow the execution of computer programs associated with the method. 
The computer programs may be implemented in software, hardware, or a combination of software and 
15 hardware. 

Further modifications and alternative enixxhrnents of various aspects of the invention will be 
apparent to those skilled in the art in view of this description. Accordingly, this description is to be 
construed as illustrative only and is for the purpose of teaching those skilled 

As such, a method for diagnosing and prognosticating the state of a system is described. Other 
20 aspects, advantages and novel features of the present invention will become apparent from the detailed 
description of the invention when considered in conjunction with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic block diagram depicting the Shannon-Weaver Model for use according 
to the invention. 

Figure 2 is a schematic block diagram depicting the information path according to the invention. 
5 Figure 3 is a schematic block diagram depicting a series of information paths according to the 

invention 

Figure 4 is a block schematic diagram depicting a computation system for implementing the 
method, according to the invention 

Figure 5 is a schematic block diagram depicting a network systernfor implementation of the 
1 0 method, according to the invention. 

Figure 6 is a block flow diagram depicting an exemplary method according to the invention. 

Figure 7 is a cross sectional view of squirrel cage induction motor. 

Figure 8 depicts Ghosh and Bhadra's [5] bond graph of a squirrel cage induction motor. 

Figure 9 depicts the stator resistances in Figure 8 redistributed to each of the stator coils. 
1 5 Figure 10* depicts a simplified representation of the signal and modulated GY element 

Figure 1 1 depicts a squirrel cage rotor with five bars. 

Figure 12 depicts a transformation of a and (3 phase currents into rotor bar currents. 

Figure 13 depicts the bond graph structure including stator and rotor bar action. 

Figure 14 depicts the bond graph equivalence used in modeling. 

20 Figure 1 5 depicts the bond graph representing stator and rotor bar action in the magnetic circuit 

Figure 16 depicts angular velocity of rotor axis and stator currents in stator winding. 

Figure 17 depicts angular velocity of rotor axis and stator currents in stator windings, at startup. 

Figure 1 8 depicts angular velocity of rotor axis and 5-currents in each rotor bar, at startup. 

Figure 19 depicts angular velocity of rotor axis and 5-currents in each rotor bar, at startup. 

25 Figure 20 depicts angular velocity of rotor axis and 5-currents in each rotor bar, from startup to 

steady state. 

Figure 21 depicts stator currents and rotor velocity of a machine with a broken rotor bar. 
Figure 22 depicts stator current of 2nd phase and rotor velocity of a healthy machine at steady 

state. 



013270.00015:124601.01 



WO 02/27418 



PCI7US01/29935 



w ~ 7 
Figure 23 depicts stator current of 2nd phase and rotor velocity of a machine with a broken rotor 

bar at steady state. 

Figure 24 depicts the angular velocity of rotor axis and 5 currents in each rotor bar when the 3*d 
bar is broken. 

5 Figure 25 depicts a torque-time plot of healthy machine and one rotor bar-broken machine. 

Figure 26 depicts rotor velocities of healthy and shorted machines. 
Figure 27 depicts rotor torques of healthy and shorted machines. 
Figure 28 depicts rotor bar currents of shorted machine. 

Figure 29 depicts Kim and Bryant's bond graph of an induction motor with state variables. 
10 Figure 30 depicts angular position and velocity of rotor axis. 

Figure 31 depicts flux in rotor & windings ; the P winding flux is s imilar . 
Figure 32 depicts flux in stator & windings ; the P winding flux is similar. 
Figure 33 depicts rotor velocity of a machine with a broken iotor bar. 
Figure 34 depicts stator current in the Frequency domain with broken bars. 

1 5 Figure 35 depicts torque-speed characteristics of the ideal and degraded machines. 

Figure 36 depicts power spectrum of the machine response and noise- 

Figure 37 depicts noise in the si gnal of the angular velocity of the degraded machine. 

Figure 38 depicts channel capacities with a broken bar. 

Figure 39 depicts rotor velocity of ideal and shorted machines. 
20 Figure 40 depicts power spectrum of angular velocity for the shorted inachine. 

Figure 4 1 depicts spectral content of stator current of phase A; (a) Ideal machine, (b) Shorted 
machine, (c) Ideal machine of [15]. (d) Shorted rnachme of [15] 

Figure 42 depicts spectral content of stator current of phase A with two severely shorted coils. 
(R s i = R S 2 =0.0079 CI, n s2 =10) 
25 Figure 43 depicts channel capacities with one shorted cofl. 

Figure 44 depicts channel capacities with two shorted coils. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Claude Shannon formulated a mathematical theory of commurdcarion. His groundbreaking 
approach introduced a simple abstraction, the communication channel consisting of a sender (a source 
of information), a transmission medium (with noise and distortion), and a receiver (whose goal is to 
5 reconstruct the sender's messages), see figure 1. 

The transmitter injects messages from an information source into the channel. The receiver accepts 
a signal from the channel that contains the transmitted signal altered by the dynamics of the channel, 
and corrupted by noise added by the channel. 

An analogy is mads between a machine component or a system and a communications channel. 

1 0 During operation, information is sent as a signal over a communications channel from transmitter to 
receiver (See Figure 2). The signal over the channel is altered by limited dynamic bandwidth, 
nonlinearities and noise. The goal is for the receiver to extract and reproduce ihe message, despite 
distortions and noise. Design of communications systems is aided by powerful theorems of Shannon 
(1949), which establish minimum signal to noise ratios for error free transmission. 

1 5 A machine component (or system) accepts a "signal" from an upstream component, by its function 

alters mat signal, and then passes the "signal" on to the next downstream component In the analogy of 
this article, a machine is a corrmiunications channel. When operating properly, the "signal" from an 
upstream component is "received" by a downstream component Faults in the machine that disrupt 
functionality alter the "signal". Faults will be viewed as agents that alter system parameters or 

20 contaminate the signal with "noise". Unless the signal to noise ratio is kept sufficiently high, 
downstream components cannot "resolve" the "signal message" error free, and the machin e 
malfunctions. 

In performing a function, a macHtc,. component, or. system accepts a stimulus "signal" from 
another upstream component, alters that signal via its mechanical- furicrion, and then passes the signal 
25 on to the next component The signal contains information, which can be envisioned as a "message" to 
other components in the machine. The "message" relates to the function or intended operation of the 
machine or machine compuiients. The mechanical function often includes kinematics of morion and 
dynamics of operation. 

Here we will strike an analogy between a machine component and a communications channel. The 
30 transmitter, an upstream component, activates our machine component "channel" with input signal x(t). 
> Vssage of the "signal" through the channel is associated with component functionality: component 
kinematics and dynamics alter the signal. The component response defines the output y(t). When the 
component channel operates properly, the "message" contained in the signal "received" by downstream 
components can be unambiguously "resolved". Defects and degradation of the component afflict normal 
35 operation, "aborting" the signal and contaminating it with "noise" n(t). Unlike electromagnetic 
communication channels, the signal may pass through multiple power domains: electrical, mechanical, 
solid, fluid, chemical, biological, etc. along its path through a component or a machine system. Wc can 
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view a machine as a channel consisting of individual component channels connected together to form a 
larger channel. 

The theory depicted and described in this application may be adapted in unique ways not 
contemplated by others, including, but not limited to, Shannon's theory. A component is designed to 
5 have functionality, which can be defined in terms of the (designer's) intended reaction of the component 
to an excitation. Degradation alters the component response. For example, fretting corrosion of the 
surface of an electrical contact changes the electrical impedance through the contact Although this 
alters the response to a voltage stimulus, the resulting signal distortions caused by changes of the 
electric contact impedance are often posed in terms of an effective noise riding on the transmitted 
10 signal. Thus degradation of the contact via fretting is often modeled as an effective noise source and/or 
an impedance change. 

In conmiunications theory, Shannon's theorems traditionally estimate the maximum rate of 
information C that can be transmitted through a coinnnmications channel, given its bandwidth w and 
ratio of signal to noise powers S/N. Designers of traditional cormmmications channels considered C to 

15 be fixed, and their designerTorts focused on designing transmission or encoding schemes that would 
increase the rate of information R up to its upper limit, the channel capacity C. If applied in a 
nontraditional manner to machinery, Shannon's theorems can yield a threshold signal to noise ratio 
(S/N)t- I* 1 the cornmunications channel analogy, dynamics inherent in the component functionality can 
be included in bandwidth w and in the signal to noise ratio S/N of the channel capacity C. These 

20 dynamics may change as the component degrades, causing C to change. For a typical machine and 
components, signal transmission rateR should be constant, since rnachine or component operation is 
often repetitive (or periodic) and at or near steady state: the machine controller and/or upstream 
components continue to inject their signals into a machine (or component), regardless of its condition. 
When applying cornmunications theory to a machine component or system, we will first trace the 

25 path of signal power flowing through a healthy (functional) component or machine, to define the 
communications channel through the component or machine. Along the signal path we will list the 
various forms of energy oi power into which the signal is transformed. Functionality will be defined in 
terms of the input to output response, for components or the system: if for a given set of input 
excitations, the output response matches within some tolerance the desired output, the component or the 

30 system is functional; otherwise it is dysfunctional. If needed, we will consider each separate 
energy/power domain and its transduction as a conmiunications channel, and then connect these 
channels together in a manner consistent with the machine's functionality and design. Bond graphs 
(Kamopp, Margolis, and Rosenberg, 1990), which map power flows through dynamic systems, can be 
useful, since bond graphs readily handle systems with diverse energy domains in an energetically 

35 consistent manner. 

After analyzing the healthy system, we will then incorporate component faults and degradation 
modes into the system model. To affect functionality, the degradation effects must alter or block the 
flow of signal power through the component. Questions we must answer include: How does each 
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degradation mode alter the signal flow, and affect system or component parameters? Does the particular 

degradation cause components to become nonlinear? Does the particular degradation generate another 

signal, Le., noise? We wffl incorporate degradation into the system model as changes to existing system 

model parameters or as additional elements (e.g., sources of noise). location of each degradation mode 

5 in the system model wffl be consistent with the locus of the degradation in the physical machine 

component or system 

Aspects of this method include: 1) Individual components, or an entire machine system consisting 
of multiple components, can be analyzed 2) System malfunctions can be predicted, mcluchng individual 
faults and those due to a collection of seemingly healthy components. Errors from slightly degraded but 
10 individually healthy components can stack through a machine system, rendering it unable to meet 
tolerance. 3) The current status of the system, and time to system maliunction can be estimated by 
simulations based on these models. 

Other aspects of the diagnostic procedure include: Determine and trace the path of the signal flow 
through the healthy system, from signal in to.signal out. For the sick system, model the faults with noise 
15 sources or parameter changes. Multiple system outputs may exist At each output, tally the signal power 
and the total noise power to obtain a signal to noise ratio S/N. Estimate the bandwidth w for the signal 
path through the degraded component communications channel, using the enhanced system modeL 
Apply Shannon's theorems to diagnose the absolute health of the machine component or machine 
system. The health of each individual component in a machine system can be assessed, and likewise the 
20 health of the entire machine system, 

The analysis of each machine component communications channel may contain the following: 
. Healthy Machine ModaL which has no faults and functions perfectly. This is an idealization that 
reflects the machine designer's original intended concept The output y 0 (t) of a signal x(t) 
propagating through this ideal machine will define the intended message or signal y 0 (t) that the 
25 machine component or machine system communications channel is supposed to transmit and 

receive. The signal powers S - P{y}and Si - P {y 0 } defined by S - P{y}= 
both in the channel capacity C and the rate of information transmission R through the channel, will 
be based on this ideal machine. The resulting model is simple, the concept of perfect health is well 
defined, and the signal y 0 (t) that the receiver is supposed to receive is well defined. 

30 

• Machine Eaulta: These include common degradation feults for a given component Common 
examples include pitting of gear teeth, fatigue cracking of shafts, and deterioration of insulation on 
electric motor stator or rotor coils. 

35 . Machine Eailtt Models incorporate the Machine Faults as sources of noise n(t) and/or changes 
in system parameters consistent with impenections, faults and degradation modes of a particular 
imchine element Noise will be defined as any signal component that should not be in the perfectly 
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transmitted and received message signal y 0 (t). This may includes harmonics generated by nonlinear 
elements. 



De graded Marhing Model : This is the overall system model that results from adding the 
5 Machine Fault Models to the Healthy Machine Model. It includes sources of noise n(t) and changes 

in system parameters. When all noise sources are zero, the healthy machine results. Transmission 
of the signal x(t) through the degraded machine (noisy cornmunications channel) induces received 
signal y*(t), generated by signal x(t) (sent through as y 0 (t))and noise n(t). 



10 The analogy may also be extended to a set of machines, a process, a rnanufacturing or assembly 

method, or others. The analogy may hold for a series of "information channels" as seen in figure 3. 

The model-based diagnostics is based on fundamental first principles of physics and 
information theory. The methods uses sensor signals to tune the parameters of a model of the system, 
such that the model then mimics the operation of the real system. Diagnostics are performed on the 

15 model. The diagnostic system can be designed as part of the design of a new machine. Also, models 

allow what if predictive scenarios for a machine's possible future health and functional condition,' given 
certain observed trends in the machine's health. For a different machine or a new design, only the 
model of the operation of the system must be altered, not the basic diagnostic algorithm. Models also 
avoids interpretation of complex sensor signals, trying to figure out what a particular peak or dip, or a 

20 band of frequencies in a signal means, in terms of machine health. Instead, time wise changes to 
machine parameters can be followed, and projection of these trends can be used to forecast future 
health. Models also permit incorporation of taiowledge of faults, and the intent of the designers of 
nracbinery, into the diagnostics routines. 

To quantitatively analyze transmission through the channel, Shannon introduced a measure of 

25 the amount of mformarion in a message. The measure is related to the probability of occurrence of the 
events for which the messages are about A message that informs the receiver that a rarely occurring 
event is about iu Lappen contains the most information. A message htforming about an already t, known" 
event conveys little information. Information entropy, a measure of the average amount of information 
(or uncertainty) in a message, can be defined [1] as 



30 # = 



-^P, log 2 Pi tf xis discrete 

(i) 

- ^ d p(x)log 1 p(x)dx if xis continuous 



Here Pi is the probability of occurrence of the message's event xj if the random variable is discrete, and 
p(x) is the probability density function for the random variable X , if the random variable jc is 
continuous. Here Pi is the probability of occurrence of the message's event*, xi if the random variable 
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is discrete, and />(*)is the probability density function for the random variable * , if the random 

variable x is continuous- 
Shannon's entropy rate (R) measured a source's information production rate, and the channel ■ 
capacity (Q measured the information carrying capacity of the channel. As per one of Shannon's 
5 theorems [1], if J* ^ C , then there exists a coding technique which enables transmission over the 
channel with an arbitrarily small frequency of errors. This restriction holds even with bounds me noise 
in the channel. A converse to this theorem states that it is not possible to transmit messages without 
errors if & > C . Thus the channel capacity is defined as the maximum rate of reliable information 
transmission through the channel. 

In another theorem, Shannon derived the channel capacity for a time continuous channel with 
additive white Gaussian noise. His expression 



10 



s 



C = <ylog 2 

involves the average transmitter power, 
1 o 

15 of the signals x o (') , the power of the noise, 
1 o 

and bandwidth CO of the channel in hertz. If the bandwidth is non-flat, then the capacity of the channel 
is given by 

" )df. (5) 



20 Similarly, the entropy or information rate for messages 

- R = co i log 2 (S i /N i ) (6) 
derived by S hanno n involves S ; , the average power of the desired signal to be transmitted, , the 
maximum allowed RMS error between recovered and original messages, and CO i , the signal bandwidth. 
Shannon's commumcation theory could be applied to tne fault diagnosis of machine systems. 

25 A machine component (or system) accepts a signal from an upstream component, by its function alters 
that signal, and then passes the signal on to the next downstream component In Bryant's analogy, a 
machine conveys information in a signal and is thus a communications channel. When operating 
properly, the signal passes through the system and is successfully received within desired tolerances at 
the machine's output. Faults that disrupt operation alter the flow of signaL Faults will be viewed as 

30 agents that contaminate the machine's signal with "noise". Unless the signal to noise ratio 
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kept sufficiently high, downstream components cannot resolve the signal message error free, and the 
machine malfunctions. 

Noise is defined as an * 'unwanted signal tending to obscure or interfere with a desired signal", 
as "any signal which interferes with the transmission of a signal through a network or tends to mask the 
5 desired signal at the output terminals of the network", and as "an unwanted signal tending to interfere 
with a required signal". Thus noise is the difference between the actual signal received, and the signal 
desired to be received. To apply this definition to mechanical systems, we must define the desired 
signal. We shall call this desired signal the "ideal" signal X 0 (f) . Note that x o (?) f an idealization, must 

be produced by a system without noise. This is possible only with models, not with real systems. 
10 The "ideal" and "degraded 11 models may be defined as follows: 

• The ideal machine model has no faults and functions perfectly. Its output defines the signal x o (0 m at 

the machine channel is supposed to receive. From this, we can estimate signal power S ( . 

• The degraded machine model is the overall system model mat results from adding faults to the model. 
We will incorporate faults as noise , Thus the signal x(t) = x o (t) + n(t) contains noise MO, 

15 defined as any signal component that should not be in the perfectly received message signal. Noise is 
any deviation from the ideal signal, including unwanted harmonics generated by nonlinear elements. 
This will estimate the noise power tf . 

We can incorporate degradation or imperfections into the system model. Degradation can be instilled 
in a bond graph model by varying bond graph parameters, adding noise (effort or flow) sources, or 
20 changing the power pathways. 

The models may take various forms. These forms may be any form appropriate for use in the system 
of application. For example, these forms may be heuristic, neural networks,-detenxiinistic, probabilistic, 
and others. 

The method and model system maybe implemented on a computer system, S see Figure 4). 

25 The term "computer system" as used herein generally describes the hardware and software components 
that in combination allow the execution of computer programs. The computer programs maybe 
implemented in software, hardware, or a combination of software and hardware. A computer system's 
hardware generally includes a processor, memory media, and input/output (I/O) devices. As used 
herein, the term "processor" generally describes the logic circuitry that responds to and processes the 

30 basic instructions that operate a computer system. The term ''memory medium" includes an installation 
medium, e.g., a CD-ROM, floppy disks; a volatile computer system memory such as DRAM, SRAM, 
EDO RAM, Rambus RAM, etc.; or a non-volatile memory such as optical storage or a magnetic 
medium, e.g., a hard drive. The term "memory " is used synonymously with "memory medium" herein. 
The memory medium may comprise other types of memory or combinations thereof. In addition, the 

35 memory medium may be located in a first computer in which the programs are executed, or may be 
located in a second computer that connects to the first computer over a network. In the latter instance, 
the second computer provides the program instructions to the first computer for execution. In addition, 
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the computer system may take various forms, including a personal computer system, m ^r r ame 

computer system, workstation, network appliance, Internet appliance, personal digital assistant (PDA), 
television system or other device. In general, the term "computer system" can be broadly defined to 
encompass any device having a processor that executes instructions from a memory medium. 
5 The memory medium preferably stores a software program or programs for the reception, 

storage, analysis, and transmittal of information produced by an Analyte Detection Device (ADD). The 
software program(s) may be implemented in any of various ways, including procedure-based 
techniques, component-based techniques, and/or object-oriented techniques, among others. For 
example, the software program may be implemented using ActiveX controls, C++ objects, 7avaBeans 5 

10 Microsoft Foundation Classes (MFC), or other technologies or methodologies, as desired A CPU, such 
as the host CPU, for executing code and data from the memory medium includes a means for creating 
and executing the software program or programs according to the methods, flowcharts, and/or block 
diagrams described below. 

A computer system's software generally includes at least one operating system such Windows 

1 5 NT available from Microsoft Corporation, a specialized software program that manages and provides 
services to other software programs on the computer system Software may also include one or more 
programs to perform various tasks on the computer system and various farms of data to be used by the 
ope rating system or other programs on the computer system The data may include but is not limited to 
databases, text files, and graphics files. A computer system's software generally is stored in non-volatile 

20 memory or on an installation medium A program may be copied into a volatile memory when rurrning 
on the computer system. Data may be read into volatile memory as the data is required by a program. 

Further, the method may be implemented across a set of networked devices (See Figure 5). 
The method may be performed remotely from the system Further, the results of the method may be 
transmitted, stored, processed, and accessed across a network, among others. 

25 For example, parameters for a model of a patients health may be stored on a smart card. 

These may be accessed ar.H.rnmbined with the method to determine a change in state of the patient's 
health. In another exemplary embodiment, a machine may be located in a remote location. A service 
provider may periodically access data from the machine from a remote location and diagnose the 
machine. These diagnoses may be used in predicting the faHure of me rnachine. Further, these 

30 diagnoses may be used in placing an order for a replacement. 

Figure 6 depicts a flowchart for diagnosing according to the invention. The method may be 
implemented in software and/or hardware. Further the method may include some or all of the steps in 
various combinations. 

In a first step, the user is directed to assemble detailed dynamic systems models of the machine 
35 system in question. The models may possess a one to one correspondence between physical components 
or elements in the real physical system, and elements in the dynamic systems model. One may inchide 
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all possible faults and potential failures in the system models. This invention may tabulate the effects of . 
faults as "noise" in the system. Noise in a signal is the difference between the actual signal and the 
expected signal. In the model, noise may be induced by changes in parameters of dynamic system 
elements, which then alters any signals passing through a system Or, if a certain fault cannot be 
5 described by these means, then sources of noise (often white noise) will be inserted into the system 
model, at locations in the model that are consistent with the locations of the faults in the real machine. 
The intensity of these noise sources can then be adjusted to make the model behave like the real 
machine. 

One may then judiciously monitor the machine or system Excite the machine or system, and 

1 0 observe the machine's or system's response, for example, via the sensor outputs. 

One may then tune the model 1 s parameters, so the model rnimics the real system Excite the 

system model with the same excitation as the previous list item Compare the outputs of the model to 

* 

the corresponding outputs of the real machine or system If the model's outputs differ from the real 
machine's or system's outputs, adjust or change values of model parameters, including the intensity of 
15 the noise sources, until the model's outputs closely match the actual system's outputs. 

One may then manipulate the model, which now mimics the real machine or system in its 
present condition: 

From the model, one may calculate the channel capacity, C, of the machine. C is the Tttaximnm 
amount of information that can be observed successfully conveyed through the machine. The channel 
20 capacity depends on the design and construction of the system, and the present condition of the system, 
which results from manufacture, aging and damage. Faults are encoded as "noise" in the model. 
Analytically, C depends on the strength of the noise levels in the system, relative to me strength of the 
excitation system response signal. 

For a desired job to be performed by the machine, one may calculate the rate of information R 
25 associated with the job R depends of the desired speed at which the machine does the job, the desired 
loads, the complexity of the job, and the desired accuracy at which the machine should do the job. R is 
measured in bits of information per second. 

- Compare R to C. If R < C, the machine will perform the desired job within the desired 
precision. If not, the system has functionally railed. 

30 The comparison of R to C may yield a diagnosis. Alternately parameters of the tuned model 

may yield a diagnosis. Further, this diagnosis may be associated with the deterrnined noise. In 
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addition, the noise and/or diagnosis maybe indicative of combined faults. Further, combined variances 

in parts, while within tolerance limits, may comprise a fault, defect, or others. 

The method may be repeated over time to build a prognosis of the machine or system For 
example, a prognosis may predict the failure of a part 
5 Further, the method may be applied to many systems such as those depicted above. In 

addition, parameters from the tuned model may indicate the type or state of a defect, fault illness, or 
condition, among others. 

In typical applications, the method may involve formulation of extremely detailed models of 
machine devices to describe a machine's condition. These are critical to success. For example, included 

10 are bond graph based models of a motor, a gear box, and other mechanical transmission components. 
These extremely detailed models: (a) exhibit a one to one correspondence between elements in the 
model and components in the real system; (b) incorporate all known effects of the device into the 
model, including defects; (c) include in the models via finite element concepts instilled into bond 
graphs the dynamically distributed nature of components in the real system, and (d) use noise sources to 

1 5 account for defects and degradation of components. Simulation of the motor and gear box models can 
generate the complex spectra measured during operation of these devices. 

The models tabulate the effects of system faults (machine maladies) in a very novel way: as 
'noise" in the machine. Noise is the difference between the actual signal received, and the expected 
signal that should be received. As a machine degrades or ages, the difference between actual and 

20 expected signals becomes larger, and thus noise levels increase. These noise methods permit 

incorporation of faults into the models that heretofore could not be described analytically. The concept 
of noise has been used heavily in electronics and communications engineering, to design around noise 
' " tr iaults" always present in these electronic and electromagnetic systems. Electronic noise, including 
resistor noise, shot noise, burst noise, and flicker noise among others has been generally tabulated or 

25 modeled with noise sources placed in a model of the electronic circuit This work imported this body of 
knowledge to mechanical, hydraulic, and other physical systems, but in addition, systems extended the 
modeling schemes of noise to include noise induced by changes in parameters of the system. 

The method also applies techniques of information theory to machinery - as opposed to present 
applications that are limited to electronic communications systems - to quantitatively assess the current 

30 health state of a machine. The method treats a machine, such as a CNC engine lathe, as a noisy 

communications channel, to assess reliability and functional condition. A message transmitted and 
received over a communications channel picks up noise due to imperfections present in the physical 
channel. For example, music transmitted over an AM channel is overwhelmed by buzzing when the 
receiver is near electrical power transmission lines: the transmitted musical message is obscured at the 

35 receiver by electrical noise. In an analogous manner, a machine transmits a message over a machine 
channel. For example, a lathe, viewed as a communications system, has transmitter = CNC controller, 
channel = (drive motor + gear box + lead screw, + tool carriage on ways + cutting tool / workpiece), 
and receiver = workpiece. "Noise" includes effects of fatigue, spurious vibration (from other machines), 
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and other errors due to wear of machine and cutting components. A transmitted ''message" is properly 
"received" if the finished part is within tolerance, or in a general machine, if the machine performs its 
function within specified tolerances. Excessive noise in the rnachine system may cause a part to be out 
of tolerance, or causes the general machine to operate outside the specified tolerance limits. With this 
5 view, Shannon's cornmunications theorems may be applied to machinery. Shannon's theorems may 
accurately estimate the limits on the amount of uifdrmation per unit time C that can be sent through a 
noisy communications channel. C depends on the channel's state, including dynamics and signal to 
noise strengths (ratios). For a lathe, making a part of certain geometric complexity at a given speed, to 
within a desired (fidelity) tolerance is characterized by an information rate R. TfTt < C, Shannon's 
10 theorems predict success; if R > C, the part will be out of tolerance. As a machine deteriorates, C 
decreases, and eventually R > C. Now the machine cannot make the part with the same speed and 
tolerance. The channel concept appears to be a very sensitive discriminator of a machine system, even 
for the stacked effects of a collection of moderately degraded components. 

The method may be used for predicting the future conditions of rnachinery, for scheduling 
1 5 maintenance and avoiding functional and catastrophic failures of said machinery. The method can 

forecast if a complex system is capable of doing a given task, within a specified tolerance. A multitude 
of parameters associated with the machine's model may be tuned, such mat the model emulates the real 
system. 

These modeling and system assessment techniques could be useful to designers of machinery, 

20 to assess the efficacy, reliability and durability of a design under various user conditions. 

In addition to mechanical systems, these methods could apply to almost any kind of dynamic 
system, including chemical, electrical, medical, manufacturing and processing, and business operations, 
among others. For example, in the medical world, a detailed model describing the dynamics of the 
caadio- vascular system could be developed. This model would possess multiple parameters that 

25 describe behavior and condition of the heart and blood vessels, and their interactions with other body 
systems such as lungs and kidneys. The model could be tuned from medical signals and data derived 
from tests and procedures, such as Electro-Cardiogram, blood pressure, and data from lau tests and 
radiology. After tuning the models, a channel capacity C could be estimated to assess the condition of 
that system, and compared to a rate of information R This comparison would assess the health state of 

30 the patient The rate of information would describe the ability of the cardio-vascular system to perform 
at various levels characterized by tasK speed, load, complexity, and tolerance. Since the rate contains 
these factors, degrees of health and sickness could be assessed quantitatively or assessed, in a formal 
manner. This could automate medical diagnostics. Medical prognostics would extrapolate trends of 
parameters in the model, or trends contained in the data, and apply the channel capacity and rate of 

35 information concepts of cornmimi cations theory, to forecast future health scenarios. 

These methodologies could be extended to evaluate business practices, procedures, and 
enterprise structures. A business operation has dynamics imposed by its processes, people, and 
structure. The application would treat an enterprise as an imperfect communications channel, and 
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construct models of information flow through that system. Transmitters— the orders— will send 
infonnation over imperfect "enterprise communication channels*. Imperfections— problems in the 
enterprise, or interference between conflicting missions-adds "noise" to channels. Receivers— the 
customers— must receive the message— the product— within tolerances— customer expectations— 
5 despite noise. The application would define "channels" through enterprise units, construct models that 
mimic these channels, and then apply conmmnications theory to diagnose and prognose these channels. 

The models in these embodiments and claims can take various forms: from structured methods 
such as bond graphs, differential equations, and finite elements, among others, to heuristic methods 
such as neural networks, fuzzy logic, expert systems, and other computer methods. 

10 The method for applying commumcation theory to machines and systems need not be limited 

to signals derived from models. The method could be extended to signals measured from real systems. 
Here the ideal signal Xo(t) could be approximated from measurements taken from a real machine, or 
from several machines, in excellent condition. The difference between Xo(t) and the signal x(t) measured 
from a degraded machine could replace those derived from models, mentioned earlier. Similarly, the 

15 difference could be used to confirm that a machine operates within tolerances. Further, an ideal signal 
could be a signal from a machine with a known defect. The difference between the signals would then 
confirm a specific defect, among others. 

F xpmpl Qr Y application ic a squirrel ra q* induction motor 

20 Equation numbers in this example refer to equations listed in this subsection. Similarly, an 

appendix is attached mat is referenced in this subsection. 

One exemplary application of the invention is a method for diagnosing an induction motor. 
- For example, a motor has two major sub systems: a rotating rotor and a sxatic stator. Induction machines 
can have a wound rotor, or a squirrel cage rotor. Widely used squirrel cage induction machines exhibit 

25 great utility for variable speed systems and are simple, rugged, and inexpensive. The squirrel cage rotor 
is a structure of steel core laminations mounted on a shaft, with solid bars of conducting material in the 
rotor slots, end rings, and usually a fen. In large machines, the rotor bars may be of copper alloy, dnven 
into the slots and brazed to the end rings. Rotors of up to 50 cm diameter usually have die-cast 
al umin um bars. The core laminations for such rotors are stacked in a mold, which is then filled with 

30 molten alurmnum In this single economical process, the rotor bars, end rings and cooling fen blades are 
cast at the same time. 

Figure 7 is a scheiratic of a squirrel cage induction motor. A substantial literature modeling 
induction motors employs Park's (1929) two-reaction theory, which accounts for magneto-mechanical 
energy transduction via multi-port inductances. From Park's model, Ghosh and Bhadra (1993) 
35 formulated the bond graph in Figure 8. We altered Ghosh and Bhadra's bond graph to partition and 
make explicit the electrical, magnetic, and mechanical energy domains; to form a one to one 
correspondence between physical components in the machine, and elements in the bond graph; and to 
append additional elements to the bond graph to make it more consistent with real induction motors. 
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When energized by an AC supply voltage, the stator coils form a radial magnetic field vector 
that rotates within the interior of the stator, about its central axis. Within this interior the stator field 
cuts through the squirrel cage rotor, including conductor bars mat extend axially. This time varying 
field induces a voltage over the rotor bars. Resulting bar currents flow hi the sequence: bar -> end ring 
_^ opposite side bar -» opposite end ring -> original bar. Induced by mis time varying current loop is a 
secondary magnetic field, which attempts to align with the stator field. However, because the rotating 
stator field induced the secondary field of the rotor, the stator field leads the rotor field, and 
consequently, the rotor chases the stator field, always following. This is motor action (Lawrie, 1987). 
The induction motor speed depends on the speed of the rotating stator field. 

The real system we will consider is a two pole, 'Y' connected three phase squirrel cage 
induction motor. In (Ghosh and Bhadra, 1993; Sahm, 1979; and Hancock, 1974), a mulri phase 
induction motor was modeled with an equivalent two axis representation. Each phase winding generates 
its own magnetic field, which can be represented as a vector aligned along the axis of the winding. The 
sum of these phase vectors produces a phasor vector. If the phase vectors vary properly with time, the 
phasor rotates. 

A transformation from three phases (a,b,c) to two phases (cc,p) was represented hi (Hancock, 
19 7 4) in matrix form. If the 4 a J and V phase windings are co-axial, the induced Magneto Motive 
Forces (MMF) of the V and 'a' phases of the three and two phase systems are co-directionaL By 
appropriate changes to the two phase currents, the magnitude of the phasors of the three and two phase 
systems can be made equal Ghosh and Bhadra (1993) represented this in their bond graph via 
transformer elements in the stator section. The two phase currents were represented in terms of three 
phases as 
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Under assumptions of a spatially sinusoidal distribution of MMFs, and ignoring magnetic losses and 
saturation, Ghosh and Bhadra (1993) expressed a symmetric indflefioh motor in an orthogonal 
stationary reference frame with a and p phases fixed on the stator as 
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(2) 



Equation (2) relates stator voltages to stator and rotor currents. In addition, needed is the electro- 
magnetic motor torque for a P-pole machine, expressed as 



30 



(3) 



This motor torque is balanced against other torques via 
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T e =J^ + cc> ln +T L (4) 

Terms on the right side of equation (4) represent rotor inertial torque, shafVbearing damping 
torque, and load torque, respectively. In equations (2) to (4), and V p> are a and p axis stator 
voltages; ias and ip s are a and 0 axis stator currents; w and i pr arc a and (3 axis rotor currents; R 5 and R r 
5 are stator and rotor resistances; I*, Lm and Lr are stator self inductance, mutual inductance and rotor self 
inductance; T c and Tl are electro-magnetic torque and mechanical load torque; J is the moment of 
inertia of the rotor, c is the viscous resistance coefficient; ©r and com are electrical and mechanical 
angular velocities of the rotor; and P is number of pole pain. 

Ghosh and Bhadra (1993) represented equations (I) to (4) in their bond graph, reproduced in 
1 0 Figure 8. They used modulated gyrators MGYrr^Wpj, MGYi 2 =Lripr, MGYrrs^Lmi^, MGYT 4 =LriaJ> 
to represent the electro-magnetic torque of equation (3); employed transformers TF:mi, TFrmz, TF:nc, 

TF:m4, 7T:m 5 with modnli^ = ^> m z = ™ 3 = -V6,ra 4 = *j2,m 5 = -4l to implement the 

mathematical transform of equation (1); and excited the system with effort sources MS e :v a , MS e :Vb ; and 
MSc:vc having sinusoidal voltages with equal amplitudes but 0, «/3, and 2tJ3 phase lags, respectively. 
1 5 Although this correctly programs the governing equations for a three phase induction motor, it lacks a 
correspondence between bond graph elements and real system components. Moreover, elements and 
their constitutive laws involve only electrical and mechanical energy dornains. Faults or design 
parameters relevant to the magnetic domain are only implicit in the mutual inductances, posed as 2 port 
inertances I:cc and 1:0 with constitutive laws 
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In Xc», \p t Xo,y and fyi are flux linkage of the respective windings. In Figure 8, five integral 
(independent) causalities exist on inertance energy storage elements, with system state variables 
/Lps, Xor, Xp, and h, where h is the rotor angular momentum 

To represent real system elements or components explicitly, certain bond graph elements 

25 should be moved, altered or added. In Figure 8, a and £ phase stator resistance elements, Rs* and R sp 
should be split into three stator coil resistances R«, R^ and R SCJ without alterating the governing 
equations. The revised bond graph shown in Figure 9 moved Rea and R^ back through the transformers 
in front of the phases. To maintain an equivalence between Figure 8 and Figure 9, we must relate R„> 
R,b and R tt to R* a and R^. Since most motors possess symmetry between phases, let R„ = Ra m Rsc = R, 

30 and R^ = Rsp = R s . For the bond graphs of Figure 8 and Figure 9 to be equivalent, the voltages (efforts) 
to the 2-port inertances on the stator sides must be equal in both Figure 8 and Figure 9. The causality in 
both Figure 8 and Figure 9 asserts that these voltages to the 2-port inertances arise from the neighboring 
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Inunctions. Summing voltages from other bonds to these 1 -junctions, and equating these respective 
voltages between Figure 8 and Figure 9 gives 



v * Ri -± 

m 4 m s m A 



m 2 772 3 



1 k-id^ 



m 2 
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(6) 



(7) 



5 By solving for i 0s /ias 3 we obtain equations in terms of resistances and transformer moduli 



(8) 



i a m 2 m i mlmlR s - {mjnjn^ + m 2 7n 3 m A )R 

mfmlnqm 4 m,R x - (m^mfo^ + D^mjn^ + m^m 4 m 5 )R 
(mfm 2 m%m 5 + n^n^m^rn^R 

By replacing the transformer moduli, m— ms of the three phase to two phase transformation with real 
numbers, m, = m 2 = ^ = m A = V2, m< = -45 , which is given in equation (1), we find that 
Rj— R, i.e., Rsa = Rsp* 51 R» = Rjb = Rsc- 

Simplified Representation of the Signal and Modulated GY Elements 



In terms of the 2-port I field of equation (5), equation (3) can be rewritten as 



(9) 



From this relation, Figure 8 can be rearranged into the form of Figure 10. where the modulated gyrators 
MGY: m = Xpr and MGY: i i2 = Aor are modulated by the flux linkages and fy r of the 2-port 
inertances. 

The number of squirrel cage rotor bars depends on the rotor's size, and usually, tens of bars 
are in one rotor. In this study we consider the squirrel cage rotor with five bars (numbered 1 to 5) 
depicted in Figure 1 1 . Shown also is the rotor magnetic field (dashed line), with north poles (N) on top 
of the rotor, and south poles (S) beneath, and bar currents. Currents directed out of plane are denoted by 
a V, and currents flowing into the plane are denoted by a V. Each end of each rotor bar is attached to 
a solid end ring. Induced currents flow through each bar and end rings. With five bars, there exist five 
different currents (flows) in this rotor. At the instant of the rotor position shown in figure 1 1 -(a), the 
sums of the currents induced by the rotating magnetic field of the stator in bar 1, 2 and 5 must be equal 
to the sum of the currents in bar 3 and 4. Likewise, the current summation of bar 1 and 5 at the position 
of Figure 1 1 -(b) must equal the sum of currents in bars 2, 3 and 4. In Figure 1 1, the thickness of each x 
and • shows the relative current magnitude in each bar. 
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To incorporate individual rotor bars into the bond graph, the a and p phase currents and 

voltages of the rotor should be split into separate bar currents and voltages. The a, b, c and cc P axes 

are stationary with respect to the stator, but because the rotor rotates relative to these axes, bar currents 

must depend on the rotation angle 6 of the rotor. Using results in Hancock (1974), rotor bar currents 

can be related to the a, P phase currents as 

iri = m[i ^ os{0+ ^^ } + i ^ M B + m^L }] (10 ) 

In equation (10), i* represents the current in the k* rotor bar (k = 1, 2, . . . n), A«r and X 9 r are rotor 
currents from Figure 8, and magnitude modulus m depends on the total number of bars, n. For n = 5 
bars, we will have currents i r i to is. Accordingly: rotor bars can be incorporated into the bond graph of 
Figure 10 via 0-modulated transformers. 

Figure 12 shows the transformation of a and (3 phase currents into individual rotor bai 
currents, where the transformer moduli are 



mr k =mcos<8 + 



L 2(ft-l)frl 
, = msm<6+ f 



k = l,2,...,n 



with n = 5. 



(H) 



(12) 



(13) 



In Figure 12, the battery of 0-junctions on the right side completes the summation of a and p phase 
currents demanded by the right side of equation (10). The voltages that sum over the two 1 junctions 
located between the I fields and the MTFs give rise to 

K + (iwiH + (mr 5 )JL s = 0 

^ + {mr 6 )X 6 + (mr^Xy + • • • + (mrj^ = 0 
Here the flux linkage \ u X2, . A.10 associated with rotor bars are located to the right of die MTF's. To 
obtain the torque contributed by each bar, equation (10) for k = 1, 2, .. ., 5 is rewritten in matrix form as 
cos 9 sin0 

COS (0 + -y) Sin(0 + y) 

cos(0 + — ) sui(0+— ) 
= m 5 5 i.e. 

cos(0 + — ) sin(0+— ) 
5 5 

cos(0 + — ) sm(0-r — ) 

The two column vectors of the 5x2 transformation matrix A form an orthogonal set for any value of 
rotor rotation angle 6; the rank of A is 2. For the mxn (mDn) matrix A having rank n, there exists 
(Strang, 1 988) an nxm left-inverse B such that BA=I D , where I n is the identity matrix of order n. In our 
model 



irotor - Aitwo phase 



(14) 



A A = m 2 



"1 0" 

0 1 



(15) 
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and the left-inverse of A is A T if m = T' i.e., the transformer modulus m has a value which 



a rotor of nbars, /» = J~ . The proof is shown in Appendix for this 



normalizes A A, For 
From equations (14) and (15), the inverse transformation is 



subsection. 



4i 
"Ts 



COS& cos(0 + — ) cos(<9 + — ) cos(0~— ) cos(0 + — ) 

5 5 5 5 

sintf sm(0-r — ) sm(0-H— ) srn(0 + — ) sin(0 + — ) 

0 J O J 



>2 



(16) 



10 



15 



If substituted into the rotor output torque equation (9), the electromagnetic torque becomes 



(17) 



The revised bond graph in Figure 13 includes stator and rotor bar interactions based on equation (17). 
Here the moduli of the k* modulated gyrator is 



(18) 



where n = 5 for Figure 13. Finally, the electric resistances of the rotor were grouped with each rotor bar 
in £ manner similar to that of the stator resistances. 

The bond graph in Figure 13 models the interaction between stator coils and rotor bars with 2- 
port I elements — inductances—in the electrical energy domain. An inductance only describes storage of 
magnetic energy. Neglected are power losses and leakage effects -in the magnetic domain, which may be 
caused by component deterioration. To describe these interactions, we replace all I inductance elements 
with equivalent combinations of gyrators and C elements, without violating causality. Figure shows 
equivalent bond graph representations between arrrand a GY and C combination; and a TF and GY 
combination. 



20 In Figure 14, n is the gyrator modulus (the effective number of coil turns); m is the transformer 

modulus; X is the flux linkage; <j> is the magnetic flux [Wb]; M is the magneto motive force fAJ; p is 
the permeance of the magnetic circuit element [H]; ei and d are efforts; and fi and {2 are flows. In 
Figure 14 -(a), through the gyrator relations A = n<p and ni — M , Using the constitutive law of the C 
element, M = §If$, the two port I elements pertaining to the a and p phases were converted into 2 -port 

25 C elements that now represent interactions between magnetic flux and magnetomotive force of the 
stator and rotor. Figure 15 shows the new bond graph with five rotor bars and the GY - C - GY 
combination that replaced the 2-port I. The gyrators were then moved through the bond graph to new 
locations more consistent with motor components. The GY to the left of the 2-port C was moved into 
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the electrical section, where it now represents the action and number of turns of the stator coils. The GY 
leap-frogged the transformers that were based on equation (1), changing moduli of these transformers 
according to Figure 14 -(b). The GY to the right of the 2-port C skipped over a 1-junction, converting 
that 1-junction into the 0-junction shown in Figure 15. Similarly, a 0- and 1-junction to the left of the 2- 
5 port C in Figure 13 were converted to a 1- and 0-junction in Figure 15. In the bond graph of Figure 15, 
electrical energy inputs, transformation of energy from electrical domain to magnetic d omain, 
mathematical phase transformations, power interactions between stator and rotor bars in terms of 
magnetic flux and magneto motive force, and mechanical rotor output are all represented and labeled 
In Figure 15, the two sets of gyrator moduli ris and Or stand for the effective coil turns which relate 
1 0 electrical and magnetic variables of stator and rotor, respectively. 

State equations were derived from the bond graph of Figure 15 with rm = ni2 = ito = n*, R» - 
Rsb = Rsc = Rs- In terms of magnetic variables, the state equations are 



K =-^-+ — + — - R,*r l W<P--PjP.) 

n s 7tx x n J n i n s**h 

for = -WH-ZjPm + P/Par) + «V ^ ^ 
9 fir = -WKrP^fi, + P.<Pfr) + <P^—7 

15 where the magnetic state variables are stator and Totor phase fluxes <)>«<, (j>ps. ^w, rjipr, and rotor angular 
momentum h. The constitutive law of the 2-port C element is 



(20) 



In the state equations, 

1 1 



(21) 



20 where the permeances P s —~\>P m = — >JL ~^ r ~ — 2" m expressed in terms of coil turns and 

n s n s n r n r 

in duc tances of stator and rotor. Here n s is the number of effective stator coil turns, n r the number of 
effective rotor coil turns, $ the magnetic flux [Wb], M the magnetomotive force [A], P is the Permeance 
[H], and hthe angular momentum [N-m-s=kg-m 2 -sec], 
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Simulations of a squinel cage induction motor used the bond graph simulation tool, 20-SIM 
(Control Lab Products, 1 998). For integration of state equations, a Runge-Kutta 4 th order method was 
adopted. Values of the system parameters for the simulations are presented in 

Table 1, some were identical to those used by Ghosh and Bhadra. 

5 Table 1 System parameters of a two pole, three phase squirrel cage induction motor 



R«, Rsb, RscMI 


Stator coil resistance 


0.07SS 


Rrl,Rc2,-..,Rrl0[Q] 


Rotor bar resistance 


0.0408 


Va, Vfc V C [V] 


Input voltage amplitude 


230 


[Hz] 


Input voltage frequency 


60 


um 


Stator inductance 


0.0153 


um 


Rotor inductance 


0.0159 




Mutual inductance 


0.0147 


n 


Number of rotor bars 


5 


Us 


Number of effective stator coil turns 


100 


Ur 


Number of effective rotor coil turns (bar) 


1 


c [N-s/m] 


Mechanical resistance 


0.15 


J [kg.m 2 ] 


Mechanical inertia 


0.4 



Shown in Figure 16 and Figure 17 are plots of rotor angular velocity and stator currents versus time. 
The rotor velocity rises slowly to a steady state value of about 377 rad/sec; the stator currents oscillate 
at the input frequency with initial large amplitude. After about 1.5 seconds, the motor reaches steady 

10 stare: the currents in stator windings decrease to a steady value and no oscillation of rotor velocity 
exists. Figure 16 plots the rotor axis angular velocity vs time when 230V, 60Hz three phase AC 
voltages are input to the stator coils. Theoretically, when 60Hz alternating inputs are given to a two 
pole AC motor, the output velocity should be 3600RPM (377 rad/sec) and me simulation yields a 
steady state value very close to this (the difference is due to the mechanical resistance load). Figure 17 

15 expands the Figure 16 time scale to show the three stator currents with 120° phase difference, during 
motor start-up. 

Figure 18-20 shows the currents in the five rotor bars and the rotor velocity. Recall there exists 2n/5 
phase difference between currents in neighboring bars. This is clearly shown in Figure 19, which 
represents the motor starting moment. While the 60 Hz frequency of the stator currents generate a 
20 constant rotational velocity of the rotating magnetic field, the frequency of currents in the rotor bars 
decrease continuously as the rotor velocity increases. This is related to 'slip' in induction motors, the 
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normalized difference between the electrical angular velocity of the aii gap MMF established by the 
stator currents, and the electrical angular velocity of the rotor (Krause and Wasynczuk, 1989). Slip is 
defined as 

a-*'"' (22) 

5 where o 5 is the synchronous speed, or the speed of the stator currents, and ©r is the speed of the rotor. 
The magnitude and frequency of the currents and voltages of the rotor depend on the relative velocity 
between the rotating magnetic field and the rotor. In these sim u la t ions, this relative velocity rnaxirnizes 
at t - 0, where the slip is unity. As the rotor velocity increases, the relative velocity and the slip 
decrease, suggesting mat the decrease of amplitude and frequency of rotor bar currents in Figure 18 are 

10 probably due to the decrease of slip. If ch - (Dr, slip s - 0 and no current is induced in the rotor bars 
(hence no torque). However, the steady state currents of the rotor bars in Figure 18 are not zero, (even 
though there is no external load) because of the Actional load of the bearing modeled as a resistance 
R:c in Figure 1 5. If an external load is applied to the motor axis, the slip should increase and therefore 
the current and voltage in the rotor bars should also increase. Figure 20 shows the currents in the rotor 

15 bars during steady state. 

AH simulation results shown above are for a healthy motor. When rotor bars break, currents, 
velocity, and torque will deviate. Because we have a one-to-one correspondence between bond graph 
elements and rnachine components, it is possible to represent broken rotor bars by increasing the rotor 
bar resistance R r . In modern squirrel cage induction motors, bars and end-rings contact the rotor core. 

20 Due to this available current shunt, currents in a broken bar are not zero (Manolas and Tegopolous, 
1997), i.e., the resistance is not infinity. Figure 21 shows the stator currents and rotor velocity for a 
rotor with the third rotor bar broken. During the transient rise time, me rotor velocity increases, and 
exhibits oscillations. Even at steady state, there exists periodic deviations of rotor velocity. -With these 
deviations, the amplitude of the currents in the stator coils also change. These changes are more clearly 

25 presented in Figure 23 ; for comparison, a corresponding healthy machine simulation is shown in Figure 
22. Figure 24 plots the currents in each rotor bar, with bar 3 assumed broken. From Figure 24, the 
induced currents are largest in the two rotor bars nearest the broken bar. Figure 25 compares the torque 
characteristics of the healthy rnachin e and broken bar rrachine. The rotor torque osculates in the broken 
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bar machine, even at steady state. During startup, the oscillation of torque is larger in the broken bar 

machine than the healthy inachiae. 

Simulations of an induction motor with a short circuited stator coil are shown in Figure 26~28. 

In these simulations, the resistance of the shorted coil decreases, and the coil current, the magnetic 

5 fields, and the induced currents in the rotor bars also change. Figure 26 shows a difference in rise time 

of rotor velocity between the healthy rnachine and the stator coil short-circuited njachine. Figure 27 

shows the rotor torque for both healthy and shorted machines. The overall trend of the torques are 

similar, but there exists small amplitude and relatively high frequency oscillations in the short-circuit 

case. These oscillations are also seen in the rotor bar currents, Figure 28, compared with the rotor bar 

10 currents of the healthy machine, shown in Figure 16. 

A bond graph model of a squirrel cage induction motor was constructed, based on a prior bond 
graph by Ghosh and Bhadra (1993), that exhibited a one-to-one correspondence between the bond 
graph elements and real system components. Included were stator coil windings for three phases, 
mathematical transforaiations to incorporate two reaction theory, magnetic state variables to represent 

15 magnetic interactions between stator and rotor, individual rotor bars and contributions to the total rotor 
torque and velocity, and mechanical inertias and resistances. The simulations in this article had five 
rotor bars. Using this model, simulations of a healthy machine were compared to simulations of 
machines with a broken rotor bar breakage and a shorted stator coil. The degraded machine simulations 
predicted oscillations in currents and angular velocities, seen in real motors. 

20 Most induction motor designs employ three phase excitation of the stator. For a rotoi with 

more bars, the bond graph of Figure 15 can be easily altered. More rotor bars can be included in Figure 
15 by adding additional pairs of power pathways to the right of the 2-port C's, such that n power 
pathways fan out from both a and p rotor phases. For the new value of w, these power pathways must 
update equations (1 1) and (12) governing moduli mi\ for the modulated transformers MTF:mr* and 

25 equation (18) governing modulus i\ of the modulated gyrators MGY:r>, To update the 
electromechanical torque, in equation (17) we must replace the 5 in the upper index of the sum and the 
square root argument in the denominator with the new value of n. 

A Ser.nnri Exemplary Application 
30 This subsection refers to equations 1-6 m the detailed description. In addition, the remaining 

equation numbers refer to equations within this subsection. Further, an appendix is attached which is 
•referenced in this subsection. 

In a further embodiment of an induction motor, the bond graph model of a squirrel cage 
induction motor from above is adjusted. This model includes stator windings for 3 phases, two-reaction 
35 theory, magnetic interactions between stator and rotor, individual rotor bar contributions to rotor torque 
and velocity, mechanical inertias, and resistances and losses. Although this model does not include 
certain critical phenomena of the induction motor - e.g., magnetic field with rotor eccentricity or rotor 
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dynamics - this model is simple and can illustrate how to apply Shannon's communication theory to 
machine systems. 

In the system shown in Figure 29, MSe:V tt , MSe:V b and MSe:V c indicate the 3-phase 
alternating voltage applied to the motor. The resistor element RiR, models resistive losses in the stator 
windings of the motor. The gyrator GYm s models the transition from the electric to the magnetic 
domain of the power flow in the system. The modulus of the gyrator n s equals the number of turns of the 
stator coil. The battery of transformers TF:m k convert the 3-phase into a rotating phasor vector. The 
two-port capacitance elements C represent the interaction between stator and rotor fields. 

In the rotor, electric voltage is induced in the metal bars by time varying flux cutting the bar 
circuits. This represented as the battery of gyrators, which have moduli n, related to the number of turns 
of the rotor. The modulated transformers MTF:mr k relate angular position of the rotor relative to the 
flux field. The resistor elements R:R, represent resistive losses in the rotor circuits. Modulated gyrators 
MGYrric convert bar currents on the rotor bars into torque; this is the magneto-mechanical interaction. 
The moduli for these gyrators depend on fluxes in the rotor. The final transformer TFim^ is related to 
ths number of magnetic poles in the system. Power lost by bearing friction is accounted for by the 
resistance R:c. The remainder of the power drives the output shaft 
Moduli in the bond graph (Figure 29) are given as follows [6]; 
1) Moduli of three phases are : 

12 111 



V 3 



2) Constitutive laws for two port C fields arc : 







M 


]■ 







(7) 



(8) 



where, ' = Q-,/5 and 



n e n r L m 



* =- 



n r 2 L, 



L s L r -Lm" ' L s L r -L m * L S L T -L m 

L s is stator self inductance, L m is mutual inductance and^ r is rotor self inductance. The gyrator 
moduli n 5 is the number of stator coil turns, and gyrator moduli n r is the mimber of rotor coil turns. 
3) The modulated transformers MTF:mrk are : 

k = l,2,...,n 



mr. 



[2 [ n , 2(*-l>r 

: i— COS< 0 + — 

V« { n 



with n = 5 



\l* { n , 

where n is the total number of bars. 
4) The moduh of the modulated gyrators MGY:r* are : 
¥ mr k - (p^mr^ j k = 1,2,...,* 



(9) 



(10) 
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5) The modulus for transformer TFmn is : 



(11) 



m m = — : P p is number of poles 

State equations were derived from the bond graph (figure 29) with n s3 = n s2 = n s3 = «r; = □= n rlC r= ?i r . 
The results of state equations are 



5 



(12) 



fcr =-( mr i 2 ^n +™r 2 2 R r2 +mr 3 2 5 r3 + W*m +""5 2 ^ 5 )^f 



-(jnr,mr tf J? rt + mr 2 mr 7 R r2 + mr 3 mr s R ri + mr 4 mr 9 R r4 + mr 5 mr l0 R rS ) — r 



1 /z 



71 r ' 772 J 



(13) 



10 



0, = -(nr 6 2 /? rl + mn?R r2 +mr l 2 R r3 +mr i 1 R r4 +mr l0 



-{pir x mr e R n + mr 2 • mr 7 ii r2 + 7nr 3 -tht^ + wr 4 • mr 9 J? r4 + 7nr 5 -iwr 10 5 r5 )— y- 

72 r 

/ \ 1 h 
-(mr 6 -Tj + mr 7 *r 2 + /7zr 8 -r 3 +mjr 9 • r A + mr l0 • r 5 ) 

71 „ '/72 IM J 



(14) 



2 A 



n s 2 + R s\ 2 J 



1 

VV 



2 \ 



6 2 , n 2 i TO 

71, +i? j2 V m 2 



772 



1 — 



R. 



2 \ 



c 2 , r> 2 



A J J 



K m 3 



m 



1 



x s2 



n.R,-, ( 1 



m,-n,{{~ n s z +R 52 2 )" n, 2 +R, 2 '{m 2 
fr 



s j 



J 



—M m + — M A 



4 



1 



m 5 • 7z, ^ 



A = (r, • mr, + r 2 ■ mr 2 + r, • wr 3 + r 4 • 7nr„ + r s ■ mr s ) — - — 

m m n r 

+(r, mr 6 + r, mr 7 +r 3 -mr, +r 4 -mr, +r 5 -mr w ) — - — -—< 

m„ -n, J 



(15) 



(16) 



(17) 



013270.00015:124601.01 



WO 02/27418 



PCT/US01/29935 



30 

where the magnetic state variables are rotor angular position Q Q and momentum h and stator and rotor 

phase fluxes 9 a* >9 ><P <* ^ d 9 fir- 

Simulation of a squirrel cage induction motor employed MATLAB®'s Runge-Kutta 4 th order 
method with a time step St - 10~ 5 seconds. Values of the system parameters presented in Table 2 

5 were given by Kim and Bryant [6, 7]. 

Using this model simulations were performed for an ideal machine, which has no faults and 
functions perfectly according to designer's specifications, and a degraded machine. The ideal machine 
will serve as a reference of desired d\Taarnic behavior. The degraded motor will exhibit common 
degradation modes, inclu ding rotor bar breakage and stator coil shorts. We will excite the ideal and 
1 0 degraded machine models with identical test signals, record these signals, and then estimate the noise as 
the difference between degraded and ideal machine responses to the same test signal. 



Table 2 System parameters of a two pole, three phase squirrel cage induction motor. 



Parameter 


Description 


Ideal 


Degraded 


Rsa> Rsb, R*c 


Stator coil resistance [n] 


0.07S8 


0.00079 - 0.0709 


Rrl>.~ JtrlO 


Rotor bar resistance [Q] 


0.0408 


0.0412-4.0800 


Vc.Vt. v c 


Input voltage anmhtude [V] 


230 




f 


Input voltage frequency [Hz] 


60 




% 


Stator Reluctance [l/H] 


5.85x10* 




*, 


Rotor Reluctance [l/H] 


563 




*. 


Mutual Reluctance [l/H] 


5.41x10* 




n 


Number of rotor bars 


5 




n s 


Number of effective stator coil turns 


100 


1-90 


n r 


Number of effective rotor coil turns (bar) 


1 




c 


Mechanical resistance [N-s/m] 


0.15 




J 


Mechanical inertia [kg.m 2 ] 


0.1 




m n (=Pp/2) 


Number of poles □ 2 


2 





15 

Figure 30~32 show sample simulation results for a nominal or ideal motor, Le., a motor 
without faults. These simulations arose from the model of equations (7) to (17), with the "Ideal" 
parameter values of Table 2. Plotted are selected motor state variables versus time, begrmung with 
motor startup, i.e. the motor voltages were switched "on". The rotor velocity rises slowly to a steady 
20 stare value of abont 377 rad/sec, as the momentum and all other state variables reach steady state. 
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"Theoretically, when 60Hz alternating inputs are given to a two pole AC motor, the output velocity 
should be close to 3600RPM (= 377 rad/sec). Rotor and stator fluxes of CC phases are shown in figure 
3 1 and 32; the P fluxes are similar. Flux amplitudes increase to steady state, consistent with the 
angular velocity. 

5 Various faults can be developed in motors. For example, stator coil shorts cause overhearing, 

increasing core losses [8]; rotor bar breaks or cracks in the die-cast rotors cause very large electrical 
resistance [6, 7, 9]; and bent or cracked shafts make the rotation wobble [10]. 

In this article, we will focus on a broken rotor bar, and shorted stator coils. When rotor bars 
break, steady state velocity and torque of the rotor will deviate from the ideal response. With the bond 
1 0 graph shown in Figure 29, a broken bar can be incorporated into the model by increasing selected rotor 
bar resistances R*. The range of deviation is given in Table 2, last column Figure 33 (upper curve) 
shows the response of the motor with a broken bar after being switched "on". In this case, the rotor bar 
resistance was increased 10 times from its nominal value of 0.0408 □,ohrns, to R, = 0.408 O.ohrns. 
Plotted is the angular velocity versus time, from start up. Here the angular velocity increases, during a 
1 5 transient time characterized by deviations of rotor velocity. 

Figure 33 (bottom curve) shows the simulated startup (step) response for a motor with a rotor 
bar having resistance increased 100 times, to = 4.08 □.ohms. This curve shows increased and 
persistent oscillations, compared to Figure 30 for the ideal machine. 

It is well established that when rotor faults occur, rotor harmonic fluxes are produced which 
20 induce currents in the stator at frequencies of Here / is the supply 

frequency, P p is the number of poles , k =1,2,3,4. . and s is slip defined as [1 1, 12] 

S = — (18) 

In equation (18) 0> z is the synchronous speed derived from the frequency of the stator currents, and 

(O r is the angular speed of the rotor. Slip can have a value from 0 to 1. 
25 Figure 34 (a), constructed by applying a Fourier transform to the steady state portion of the 

simulation results of Figure 33 (bottom curve), shows some frequencies of the stator current of phase A 

in the vicinity of the excitation frequency (60Hz). Figure 34 (b), scanned from reference [13], shows 

spectral densities of typical currents versus frequency measured from a motor with three broken bars. 

Comparison of Figure 34 (a) and (b) show similar shape and location of spectral peaks. 
30 Figure 35 compares the torque characteristic of the ideal machine and broken bar machine R r = 

0.408 □.ohms. The rotor torque oscillates whenever the rotor velocity oscillates. Due to the rotor 

asymmetry the level of pulsating torque is increased [6, 12]. 

The average power in a signal x(t) f of duration T can be estimated as [14] 

■* 0 
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or as 

If x„ is a sequence sampled from *(*) at equally spaced discrete instants. The power spectral density, 
the magnitude squared of the Fourier transform of x(t) , is given by 

5 s(/H*0f (21) 

-tu 

For discrete x„, we employed a last Fourier transform to obtain Xu. 

The total power can be calculated in the frequency domain, or in the time domain by Parscval's 
theorem. [14] 

,o )\ X (t)\ 2 dt=)\X(ffdf (22) 

The discrete form of Parseval's theorem is defined as [14] 



Equation (5) can be rewritten as 

is c= J{iog 2 (s+;\o-iog 2 w}<to. (24) 

o 

Combination of the original signal power spectral density (S) and the noise power spectral density (Af) 
represents the signal power spectral density 

S.=S + N ( 25 ) 

from tiic degraded machine. 
20 Shannon [I] assumed a Gaussizai white noise statistically independent of the signal. To remove 

this restriction, we need to calculate the noise power directly from the time domain signals. In the time 
domain, the noise is defined as the difference between actual and ideal signals 

«(0 = x(t)-X a (t). (26) 
Here x(t) is the output of the degraded machine, and x 0 (0 is the ourput of the ideal machine. As 
25 demonstrated in the Appendix for this subsection, removal of the independence restrictions between 
signal and noise admits negative channel capacities. 

Power spectral densities K and N' can be defined as the magnitude squared of the Fourier 
transforms for signal x(t) and noise "( r ) respectively. To calculate the channel capacity with these 
values, we must replace {S + N) in equation (24) with Si , and N with N' , to have 
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c= ]{ l08 {v)} da> - (27) 

Figure 36 (upper line) shows the power spectrum of the rotor velocity of the ideal machine as 
shown in Figure 30 (upper line), and defined in section 4.1. This figure was constructed by applying a 
fast Fourier transform to the angular velocity data of Figure 30 (upper line). In this case, we assume 
5 zero noise, and thus the system functions perfectly, according to the design specifications. Using 
equation (27), we get an infinite channel capacity for an ideal machine system, since by definition, the 
noise and noise power are zero. 

If there are faults such as broken rotor bars as mentioned in section 4.2earlier, the power 
spectrum will change as noise contaminates the signal. Figure 37 shows the startup response x(t) of 
10 the machine with a broken bar (upper line), and the noise in the time domain from the degraded 
machine, defined by equation (26). This noise (lower curve and magnified in Figure 37) is the 
difference between the degraded machine's response curve in Figure 37, and the startup response of the 
ideal machine in Figure 30. The presence of several frequencies is evident. Figure 36 (dots in the upper 
line) shows the power spectra (signals and noise) of a degraded rnachine, with a cracked rotor bar. Note 
15 that the power spectrum of the degraded machine signal x(t) nearly overlaps the power spectrum of 

the ideal machine signal * 0 (0 ; in the figure, the two almost coincide. In our model, we increased the 
resistance of broken (cracked) bar by 10%. Using equation (27), we obtained a channel capacity of 
1.1 xlO 6 (bits per second). Here the integration bandwidth <*> in equation (27) was equated to the 
entire sampling bandwidth (5xl0 4 ifc) based on the Nyquisfs sampling rate, where St = 10 5 

20 seconds was the time step employed in the numerical solution routine, m this procedure, we viewed the 
numerical solution's data points as a "sampled" signal, with sampling interval equal to the numerical 
method's time step. The Nyquist's sampling rate gives the smallest bandwidth associated with the 
sampling interval St -10 5 seconds. 

In equation (6) for entropy rate R, S; represents the average power of the output signal fiom 

25 the healthy rnachine and A, represents the largest acceptable deviation, i.e., a tolerance on the noise. 
The signal bandwidth (©*) was equated to G>\ see the previous paragraph for justification. The 
functional requirements of the machine determine the noise or error tolerance N i demanded by the 

machine to work satisfactorily. For example, if we have an application wherein the maximum allowed 
error or tolerance must be within 10% of the signal of the ideal machine, and if we employ the same 
30 bandwidth as for the channel capacity, then from the equations (6) and (19), the information rate (R) is 



f 



R = a> t log 2 (S i I N,) = 5 x 1 0 4 log 2 



= 3.3xt0 5 (bits per second) (28) 
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With a channel capacity of 1.1 *10 6 (bits per second), the result for R above satisfies the condition of 
R<C . This suggests a still "healthy" machine. 

If the resistance of the broken bar increases 21 times to Rr - 0.8568D, then the channel 
capacity drops to 2.4xl0 5 (bits per second) below the required (R) of 3.3xl0 5 (bits per second). 
5 Since this result doesn't satisfy R^C, according to Shannon's theorem, the machine is 
malfunctioning. As the magnitude of the fault (bar resistance) increases, the channel capacity 
diminishes. 

Figures 38(a)-3S(d), show selected power spectral densities of the stator current of phase A at 
steady state, for selected bar resistances. These figures are similar to Figure 34 (a) and are often used as 

10 diagnostic indicators. Side bands are absent for smaller values of but start to appear after the rotor 
bar resistance equals approximately 0.7670Dohms (1780%). Figure 38(e) also plots the channel 
capacities versus the percent change of bar resistance from the ideal value given by table 2. The clashed 
line indicates the 10% noise power tolerance {R =3.3x 10 s ) estimated in the previous paragraph. Here 
percent change is defined as VoR^ = {R bar " R ° bar ) /R l !r > where R bar is the current value, and 

15 /?^r the ideal. The channel capacity of point (d) in figure 38(e) has a negative value; see the Appendix 
for reasons. As shown in figures 38(c) and 38(d), significant side bands with large intensity begin to 
appear in the power spectra, wherever the channel capacity curve sinks below the 10% information rate 
line (dashed curve). From a practical standpoint, for industrial grade rugged machin es such as motors, 
we would begin to notice errors when these exceed 10% or more in the motor's output velocity. Thus, 

20 1 0% was chosen as the critical velocity condition. 

The curve of figure 38(e) can be separated into regions with three distinct slopes: region 0, 
which connects the infinite channel capacity of the ideal system to that of "real" systems; region I, with 
stable C and "healthy " operation (region I would be associated with the normal life cycle operating 
region of the system); and region H where C declines to the (dashed) failure line. Note that the marked 

25 change in the slope of C or the rapidly diminishing values of C, going from region I to ^ could presage 
iailure. Figures 38(b)-(d) suggest that once side bands appear, the slope of C becomes noticeably more 
negative. 

Signal based diagnostic methods, that trigger upon detection of side bands, at earliest would 
notice the broken bar fault at point (b) in figure 38(e); figure 38(b) suggests that detection of the tiny 
30 side band would be difficult. In contrast, the channel capacity curve's knee - where the slope abruptly 
changes - occurs at 1500%, before 1780% of figure 38(b). Here the abrupt change in slope might be 
easier to detect 

Simulation of an induction motor with short circuited tarns on its stator coil is shown in figure 
39. Here stator resistance R.i of phase A was decreased 50%, and the effective number of turns 
35 represented by gyrator modulus n s was similarly decreased from 100 to 50. In this simulation, only one 
of the stator coils has shorted turns. In the model and physically, as turns are short circuited, the 
resistance in that coil decreases, and the effective number of turns also decreases. Figure 40 shows the 
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various signal and noise power spectra of this shorted machine. Note again that the power spectra of the 
ideal and degraded machines nearly coincide. The shorted coil seems to affect the angular velocity 
relatively less than a broken bar. The two startup responses in the figure 39 are nearly the same, 
especially at steady state. Essential differences at steady state are in the signal's phase, generally not 
5 contained in power spectra. By taking the difference between degraded and ideal response, the 
information on phase differences is conveyed in the noise power spectrum, in addition to the magnitude 
information. 

When stator coil turns short out, we observe only a rise in some of the frequency components 
which already exist in the stator ctrrrent spectra of an ideal inachine [15], Figures 41 (a) and (b) shows 

1 0 spectral content of the steady state stator current of phase A, from simulations of the bond graph modeL 
For comparison, spectra from Gojko and Penman's model [15] are also shown as figures 41 (c) and (d). 
Figure 42 shows spectral content of stator currents for two shorted coils, phases A and B. In figures 43 
and 44 are plotted the channel capacities versus percent change in the coil resistance, for shorting of 
phases A, and A and B, respectively. Again rnforrcation rate for the 10% noise level on angular velocity 

15 is shown as the dashed line 

Similar to the broken bar case of figure 38, figure 43 exhibits a 'Taealthy" region I. with stable 
channel capacity, and a region II with sharply dirninishing channel capacity. Again the sharply changed 
slope of Tegion II could prognose failure. 

In this article, it was demonstrated how Shannon's theory of coinrramication could be applied 

20 to rriachinery, to utilize Shannon's powerful theorems. Concepts of rate of information rate and channel 
capacity were reviewed, and applied to an induction motor. At the heart of the method is machine 
"noise", estimated as the difference between actual and ideal responses. By subtracting the ideal 
response x c (0 from the actual response x(t) t the noise signal contains only information about the 
faults. From the noise and signal were calculated power spectra, used in equation (27) for channel 

25 capacity. Rate R from equation (28) serves as a critical values dependent on the system's tolerance to 
errors, here called "noise". The channel capacity was calculated for a motor with shorted stator coils 
and broken bars. The channel concept agreed with other existing fault monitoring methods, but results 
suggest that it could detect faults much earlier. It can be concluded from this study that the channel 
capacity concept could serve as an effective discriminator of motor and machine system health. 

30 It is to be understood that the forms of the invention shown and described herein are to be 

taken as the presently preferred embodiments. Elements and materials maybe substituted for those 
illustrated and described herein, parts and processes may be reversed, and certain features of the 
invention may be utilized independently, all as would be apparent to one, skilled in the art after having 
the benefit of this description of the invention. Changes may be made in the elements described herein 

35 without departing from the spirit and scope of the invention as described in the following claims. 

As such, a method for diagnosing the state of a system is described, in view of the 
above detailed description of the present invention and associated drawings, other 
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modifications and variations will now become apparent to those skilled in the art. It should 

also be apparent that such other modifications and variations may be effected without 

departing from the spirit and scope of the present invention as set forth in the claims which 

follow. 
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Nomenclature 

C channel capacity 

c viscous resistance coef&cient 

/ supply frequency 

information entropy 

h angular momentum [N*m-sec] 

/ moment of inertia 

Ls, Lm, Lr statbr self inductance, mutual inductance and rotor self inductance 

M , M p a and J3 axis magneto motive force 

10 m j~m s moduli of three phases 

n modulus of gyrator (number of coil turns) 

n s number of effective stator coil turn 

n r number of effective rotor coil turn 

w(0 noise in time domain 

15 P P number of poles 

p probability of occurrence 

R entropy rate 

Rs, Rr stator and rotor resistances 

r i^rs modulated gyrator moduli of rotor 

20 51,, 51 B( sxator reluctance, mutiial rehictance, rotof reluctance 

S, N average power of the signal and noise 

5* average power of the signal including noises 

s slip 

V a , Vb, V c sinusoidal input voltages 

25 x(t) output of the degraded m^^bme in time domain 

x o (0 output of the ideally healthy machine in time domain 

CO bandwidth 

P e 9 *Pfi ot and J3 axis fluxes 

O flux linkage 
30 V^, V ps a and p axis stator voltages 

ias, ip s a and p axis stator currents 

W, ip r a and p axis rotor currents 

Rs, Rr stator and rotor resistances 

L,, 1^, L, stator self inductance, mutual inductance and rotor self inductance 

35 T t , Tl electro-magnetic torque and mechanical load torque 

J moment of inertia 
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c viscous resistance coefficient 

Q xt fflm electrical and mechanical angular velocities of the rotoi 

p number of pole pairs 

X flux linkage 

5 mi-m 5 moduli of transformers for -J phase to 2 phase transformation 

V a ,V b ,V b sinusoidal input voltages 

mi~ai5 transformer moduli 

j rk current in the rotor bar, 

m magnitude modulus that depends on the total number of bars 

1 o n modulus of gyrator (number of coil turns) 

<j> magnetic flux [Weber (Wb)] 

M magneto motive force [Ampere (A)] 

-p permeance of circuit element [Henry (H)] 

R reluctance of circuit element [1/Henry (H" 1 )] 

15 ei,e2 effort 

fufz flow 

n 3 number of effective stator coil turn 

Dl number of effective rotor coil turn 

h : angular momentum [N-m- s=kg*m 2 *sec] 

20 Rs_alpha, 
Rs_beta, 

Rsb, Rsc, R electrical resistances 



25 
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APPENDIX FOR FIRST SQUIRREL CAGE INDUCTION MOTOR EXAMPLE 



10 



Iii this section, for a rotor with n bars, we prove 



(A.1) 

5 From equations (15) and ( 1 6), the transformation matrix times its transpose is 



A A = m 



COS0 

' In 
6 + — 
n 

(* 

cos 0 + — 
n 



^A (n 4*0 (n . 2(*-l)*Q 

6 cos \6 + — cosl G+ — ... cos \0 + — — 

a . (n *A • (n 4 A ■ (* ■ 2(n-l)A 

sin (9 sin|0 + — I sin|^+ — J — sm|^0 + ~ —J 

sine? 



(A.2) 

The result of the multiplication is a square matrix of dimension 2, 



A T A = m 2 



(A-3) 



where 



t-r I 

'22 =Z sin2 l 



2(*-l)?r 



1-1 



0 + 



(A.4) 



Equations (A.4) can be rewritten using double angle trigonometric formulas: 
n^ 

15 



— -J ^ I cos{29+j7} 

W z j-o n 



312=7 S sm{2G-H- u ^-}-~ 2 sm{2e+j— } 
k=l j=0 31 

n l-co S{ 2 0 + ^ } j 
k=l j=0 



S22= I~ 



(A.5) 



2 2 2.^""*'" "n 

Via formulas 1.341-1 and 1.341-3 in Gradshtcyn and Ryzhik (1980), the sums of sine and cosine terms 
on the right sides of equations (A.5) are zero, for n > 3. Thus 
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S„ = 



*2! = 



f 1 1 

- + i-cos2<9 
2 2 
1 + cos 29 
n 

2 

(A.6) 

i-icos2<9 
2 2 
1- cos 20 

« 
2 



;n = i 

;n = 2 
;n^3 



;n = i 

;n = 2 
;nS3 



(A.7) 



5,5 =^i=1 



-Uia20 
2 

sin 20 



;n = l 
;n = 2 
;n>3 



(A.8) 

Taerefare, we can conclude 



A T A = m' 



L J 21 J 22. 



= m~ 



(A.9) 



2j 



,»ri oi 

=m l[o ij 



10 for the rotor with more than 2 rotor bars, i.e. 5 n> 3. 
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Appendix For Second Example of Induction Motor 

Range of Channel Capacity Values 

In equation (27), average power of the output signal including noise is defined as 
1 o -* c 

5 =4 Jk(0j *4 It 1 ®!**? 

o i o -* 0 

(A.1) 

Using equations (3) and (4), we get 
2 r 
* 0 

(A.2) 

10 Shannon assumed statistical independence of X 0 (0 and "(0 , which made the last term of equation 
(A.2) vanish. In this article, wc will remove this restriction, allowing for forms of noise 
n(O = -K-v„(/),(0<K<l) 

(A3) 

that can extinguish the signal such that 
15 x(O = ^(0 + «(O = (l-Kk(0 <x 0 (t). 

(A.4) 

Here the amplitude (l - K) of x(t) dhninishes with increasing noise: tins can affect the signal power 
significantly. With the values just derived, 




20 (A.5) 

can be less than unity, and the channel capacity can have negative values: as noise power proportional 
to K 2 increases, the output power proportional to (l ~ K) decreases. 
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We claim: 

1 . A method of modeling a mechanical system comprising a plurality of physical components, 
comprising: 

preparing a model of a mechanical system in which at least a portion of the physical 
5 components of the mechanical system are individually modeled, wherein the model is 

configured to output data representative of the condition of the mechanical system in response 
to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predetermined operating 
conditions during use; 

1 0 modifying the model such that the outputted data of the model in response to the 

predetermined operating conditions is representative of the condition of the mechanical system 
in response to the predetermined operating conditions. 

2. The method of claim 1; further comprising predicting a failure time of the mechanical system 
using the modified model 

15 3 . The method of claim 1 , wherein all the physical components of the mechanical system are 
individually modeled. 

4. The method of claim 1, wherein a possible fault for each of the individually modeled physical 
components are incorporated into the model. 

5 . The-kielhod of claim 1, wherein a potential failure for each of the individually modeled 
20 physical components are incorporated into the model. 

6. The method of claim 1, wherein a plurality of possible failures for the individually modeled 
physical components may interact, rendering failures not specifically associated with any single 
component, but arising from interactions between components. 

7. The method of claim 1, wherein the condition of the modeled mechanical system is 
25 represented within the model as noise. 

S . The method of claim 1, wherein the condition of the modeled mechanical system is 
represented within the model as noise, and wherein the condition of the modeled mechanical system is 
determined by calculating a signal to noise ratio for the model. 

9. The method of claim 1, further comprising: 
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calculating the channel capacity of the modeled mechanical system, wherein the channel 
capacity is representative of the design of the system and the present condition of the 
mechanical system; 

calculating a rate of information for a predeterrnined job to be performed by the mechanical 
5 system, wherein the rate of information is representative of the speed, loads, complexity and 

desired accuracy of the job; and 

s 

comparing the rate of information to the channel capacity, wherein if the rate of information is 
less than or equal to the channel capacity the model will output data indicating that the 
mechanical system is capable of performing the job at the appropriate speed, load, and 
10 accuracy. 

1 0. A computer implemented method of modeling a mechanical system comprising a plurality of 
physical components, the method comprising: 

preparing a model of a mechanical system in which at least a portion of the physical 
components of the mechanical system are individually modeled, wherein the model is 
1 5 configured to output data representative of the condition of the mechanical system in response 

to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predetermined operating 
conditions during use; 

modifying the model such that the outputted data of the model in response to the 
20 predetermined operating conditions is representative of the condition of the mechanical system 

in response to the predetermined operating conditions. 

11. A carrier medium comprising computer instructions, wherein the program instructions are 
computer-executable to implement a method of modeling a mechanical system comprising a plurality of 
physical components, the method comprising: 

25 preparing a model of a mechanical system in which at least a portion of the physical 

components of the mechanical system are individually modeled, wherein the model is 
configured to output data representative of the condition of the mechanical system in response 
to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predeterrnined operating 
3 0 conditions during use; 
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modifying the model such that the outputted data of the model in response to the 
predetermined operating conditions is representative of the condition of the mechanical system 
in response to the redetermined operating conditions. 



5 12. A method for diagnosing a state of a system the method comprising: 
measuring a signal from the system; 

comparing the signal from the system and an expected signal to determine a noise signal 
associated with the signal from the system; 

determining a signal strength associated with the signal from the system; 

1 o determining a rate of information, the rate of information associated with a desired operabihty 

of the system; 

determining a channel capacity from the noise signal and the signal strength, the channel 
capacity being a function of a frequency spectrum of the signal; 

comparing the rate of information to the channel capacity to diagnosis the state of the system, 

15 . 

13. The method of Claim 12 wherein the expected signal is a sig na l measured from an exemplary 
system operating in a known state. 

14. The method of Claim 12 wherein the expected signal is the output of a model. 

20 

15. The method of Claim 14 wherein the output of the model is adapted to approximate the measured 
signal. 



16. The method of Claim 12, the method further comprising: 
25 repeating the steps of the method over time to determine a set of diagnoses, 

deterrnining a prognosis of the system from the set of diagnoses. 
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